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Mott physics and spin fluctuations: a unifled framework 
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We present a formalism for strongly correlated electrons systems which consists in a local approxi¬ 
mation of the dynamical three-leg interaction vertex. This vertex is self-consistently computed with 
a quantum impurity model with dynamical interactions in the charge and spin channels, similar 
to dynamical mean field theory (DMFT) approaches. The electronic self-energy and the polariza¬ 
tion are both frequency and momentnm dependent. The method interpolates between the spin- 
fluctuation or GW approximations at weak coupling and the atomic limit at strong coupling. We 
apply the formalism to the Hubbard model on a two-dimensional square lattice and show that as 
interactions are increased towards the Mott insulating state, the local vertex acquires a strong fre¬ 
quency dependence, driving the system to a Mott transition, while at low enough temperatures the 
momentum-dependence of the self-energy is enhanced due to large spin fluctuations. Upon doping, 
we find a Fermi arc in the one-particle spectral function, which is one signature of the pseudo-gap 
state. 

PACS numbers: 


Strongly-correlated electronic systems like high- 
temperature cuprate superconductors are a major chal¬ 
lenge in condensed-matter physics. 

One theoretical approach to cuprates emphasizes the 
effect of long-range bosonic fluctuations on the elec¬ 
tronic fluid, for example long-range antiferromagnetic 
(AF) fluctuations due to a quantum critical point [1- 
6]. These bosonic fluctuations are also central to ap¬ 
proaches such as the two-particle self-consistent approx¬ 
imation (TPSC [7-11]), the GW approximation [12] and 
the fluctuation-exchange approximation (FLEX [Id]). 

Another approach focusses, following Anderson [14], 
on describing the Mott transition and the doped Mott 
insulator. In recent years, dynamical mean field the¬ 
ory (DMFT) [15] and its cluster extensions like CDMFT 
[16, 17] or DCA [18-20] have allowed for tremendous the¬ 
oretical progress on the Mott transition both for models 
and realistic computations of strongly correlated mate¬ 
rials [21] . In particular, numerous works have been de¬ 
voted to the one-band Hubbard model, mapping out its 
phase diagram, studying the d-wave superconducting or¬ 
der and the pseudogap [22-45]. Cluster DMFT is one 
of the few methods designed for the strong-interaction 
regime to have a simple control parameter, namely the 
size Nc of the cluster or the momentum resolution of the 
electronic self-energy. It interpolates between the DMFT 
solution {Nc = 1) and the exact solution of the Hub¬ 
bard model {Nc = oo). Despite its success, this method 
nonetheless suffers from severe limitations: i) it does not 
include the effect of long-range bosonic modes of wave¬ 
lengths larger than the cluster size; ii) the negative sign 
problem of continuous-time quantum Monte Carlo has 
so far precluded the convergence of the cluster solutions 
with respect to Nc in the most important regimes like 
the pseudogap; in) the k-resolution of the self-energy is 
still quite coarse in DCA (typically 8 or 16 patches in 
the Brillouin zone, see e.g. [31, 33, 45, 46[), or relies 


on uncontrolled a posteriori “periodization” techniques 
in CDMFT [17[. 

Several directions beyond cluster DMFT methods are 
currently under investigation to address these issues, such 
as GW + DMFT [47-53], the DTA method [54-57], the 
dual fermion [58] and dual boson methods [59, 60], or 
combinations of DMFT with functional renormalization 
group methods [61] . 

In this letter, we discuss a simple formalism that uni¬ 
fies the two points of view mentioned above. It is de¬ 
signed to encompass both Mott physics a la DMFT 
and the effect of medium and long-range bosonic modes. 
It interpolates between the atomic limit in the strong- 
interaction regime and the “fluctuation-exchange” limit 
in the weak-interaction regime. It consists in decou¬ 
pling the electron-electron interaction term by Hubbard- 
Stratonovich bosonic fields and making a local self- 
consistent approximation of the lattice’s electron-boson 
one-particle irreducible vertex, using a quantum impu¬ 
rity model similar to the one used in DMFT. It can be 
formally derived from a functional of the vertex given by 
three-particle irreducible diagrams [62, 63] . In the fol¬ 
lowing, we will therefore denote this method as a triply- 
irreducible local expansion, or TRILEX. Already at the 
single-site level, it produces, in some parameter regimes, 
a momentum-dependent self-energy and polarization, at 
a small computational cost, similar to solving Extended 
DMFT (EDMFT) [64-66]. In the following, we first in¬ 
troduce the method; we then present the solution of the 
single-site version of TRILEX for the two-dimensional 
Hubbard model. 

We focus on the Hubbard model defined by the follow¬ 
ing Hamiltonian: 

^ = X ^ X (1) 

ij i 

The indices i,j denote lattice sites, cr =t,4'i Cic 
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(the Fourier transforms of -{ciarCj^Q) and 
^/"lattice quantities'\ respectively) are given by Dyson equations: 



Figure 1: (color online) Description of the TRILEX method 


are fermionic creation and annihilation operators, and 
riia = cl^Cia- tij is the tight-binding hopping matrix 
{tij = t{t') for (next-)nearest-neighbors), while U is the 
on-site Coulomb interaction. We rewrite the operators of 
the interaction term as: 


U i ^ nf (2) 

I 

where is the bare interaction in channel I, and 
= Scrcr' where = 1 and are 

the Pauli matrices. In this paper, we consider two de¬ 
couplings: (a) in the charge and vector spin channel 
(/ = 0,x,y,z) (“xi/z-decoupling”), 11^ = 11^ = 11^ = 
and C"P and satisfy: U = - 3C/"P; 

(b) in the charge and longitudinal spin channel only 
(/ = 0, z) (“z-decoupling”), U = In both cases, 

we have two channels, denoted as 77 = ch, sp. In this pa¬ 
per, we fix the ratio to = U/2 and ?7 ®p = — [//6 
{xyz decoupling) and = U/2 and C/®p = —U/2 (z 
decoupling). We now decouple (2) using real bosonic 
Hubbard-Stratonovich fields i/f (r) in each channel and 
at each lattice site, so that the action becomes: 
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and Ciar are conjugate /3-antiperiodic Grassmann 
fields, and A'^ = 1. We are now dealing with an interact¬ 
ing lattice problem with a local electron-boson coupling. 
The lattice Green’s functions G(k, fw) and IF''(q, iff) 


G(k, iw) = [fw-I-/r — e(k) — E(k, 7w)] ^ (4a) 

IF''(q,fD) = G''[l-C/''P’'(q, 7 D)]”^ (4b) 

k and q are momentum variables, iui{i^l) stands for 
a fermionic (bosonic) Matsubara frequency, e(k) is the 
Fourier transform of t^-, and y is the chemical poten¬ 
tial. The fermionic and bosonic self-energies E and P'' 
are given by the exact expressions (written here for the 
paramagnetic normal phase) (see e.g [67]): 

E(k,7c.) = - ^ (5a) 

q,in, lUJ+lVl 

77=ch,sp 

P''(q,fD) = 2^A''G q+k. Gk.^.A^k^q, (5b) 

k,iu] iui,ih 

Here, rrich = 1, Wsp = 3 {xyz decoupling) or rn-sp = 1 
{z decoupling). A’'(q, k, zw, zH) is the exact one-particle 
irreducible electron-boson coupling (or Hedin) vertex, 
namely the effective interaction between electrons and 
bosons renormalized by electronic interactions. 

The main point of this paper consists in approximat¬ 
ing the vertex A''(q, k, zw, zH) by the local, but two- 
frequency-dependent A[(^p(za;, zH) computed from a self- 
consistent quantum impurity problem: 

A''(q, k, iu}, in) « A[[^p(iw, zH) ( 6 ) 

This strategy radically differs from DMFT, EDMFT and 
GW+DMFT which approximate the self-energy E (and 
P), not A. It implies that our E and P (computed 
from (5a-5b)) are, in some parameter regimes, strongly 
momentum-dependent while containing local vertex cor¬ 
rections which will be essential to capture Mott physics 
(see also [50]). Formally, DMFT is a local approximation 
of the two-particle irreducible Luttinger-Ward functional 
‘hLw[G'] [15, 68 [. In contrast, our approximation can be 
defined as a local approximation of the three-particle ir¬ 
reducible functional /C[G, VF, A] introduced in [62, 63] as 
a generalization of 4 >lw to higher degrees of irreducibil- 
ity. We therefore denote it as TRILEX, triply-irreducible 
local expansion. It makes it exact in the limit of infinite 
dimensions. The formal derivation of the method will be 
provided elsewhere [69] . 

The action of the impurity model reads: 

-S'imp = ~ jj dTdT'^cl^g{T - T')c^r' 

dTdT'n[u\T -T')n[, (7) 

j JJ 0 

This is an Anderson quantum impurity with retarded 
charge-charge (/ = 0) and spin-spin (J = x,y,z in 
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Figure 2: (color online) (T, U) Phase diagram (half-filling, 
t' = 0). The black diamonds (resp. squares) denote Tap 
( resp. the dashed lines are guides to the eye. The 

green diamonds (resp. squares) denote converged TRILEX 
solutions in the 2 (resp. x, y, z) decoupling. A, B, and C are 
defined as A: fUD = 96, U/D = 0.5, B: /3D = 24, U/D — 2, 
C: PD = 48, U/D = 4. The red dotted line denotes 
for the square lattice (from [70]). 


the xyz-decoupling, I = z in the z-decoupling) inter¬ 
actions. The bosonic fields (j)^ have been integrated out 
to obtain a fermionic action with retarded interactions 
amenable to numerical computations. We compute the 
fermionic three-point correlation functions to reconstruct 
the electron-boson vertex Aimp (as shown in the Suppl. 
Mat., section B). Finally, Q and derive from the self- 
consistency conditions as follows: 


Q ^(fw) = Gj^(,(fw)-h Eioc(za;) (8a) 

= TO-'(*a;) + PlM (8b) 

where, for any X, Arioc(*w) = X(k, fw). At conver¬ 
gence, this ensures that Gioc = Gimp and 
and the susceptibility are related by: 

w’'(q, in) = u^ - c/''x''(q, in)u^ (9) 

The computational scheme is illustrated in Fig. 1. 
From the impurity electron-boson vertex Aimp, we com¬ 
pute E(k, iuj) and P’’(q, in), which are then used to com¬ 
pute Q and for (7). We solve the quantum impu¬ 
rity model exactly by a continuous-time quantum Monte- 
Carlo algorithm [71] in the hybridization expansion [72] 
with retarded density-density [73] and vector spin-spin 
interactions [74]. The computation of the three-point 
functions are implemented as described in [75]. We iter¬ 
ate until convergence is reached. Our implementation is 
based on the TRIQS library [76]. 



Figure 3: (color online) Left: Evolution of the local vertex 
Re A’’(lain, iff m) (half-filling, t' = 0). A, B and C are defined 
in Fig. 2. The dashed lines denote the atomic vertex A^^ 
(Eq. (Al)). Right; ImEioc(*Wn) for TRILEX and DMFT 
(paramagnetic phase). 


TRILEX provides a unified framework for spin- 
fluctuation approaches and Mott physics. Indeed, (i) 
at small interaction strengths, the local vertex reduces 
to the bare, frequency-independent vertex X so that 
E is given by one-loop self-consistent diagrams, as in 
spin fluctuation theory in its simplest form (spin channel 
only), the GW approximation (charge channel only), or 
in FLEX limited to particle-hole diagrams; similarly, P’* 
becomes equal to the “bubble” diagram; (ii) it is exact in 
the atomic limit (t = 0): the effective local action turns 
into an atomic problem, A into the atomic vertex Aat 
(Eq.Al), and E and P become local, atomic self-energies. 

Let us now apply the TRILEX method to the Hubbard 
model on a square lattice. All energies are given in units 
of the half-bandwidth D = 4|t|. The Brillouin zone is 
discretized on a 64 x 64 momentum mesh. We restrict 
ourselves to the paramagnetic normal phase. 

In Fig. 2, we present the phase diagram in the (T, U) 
plane at half-filling. We obtain converged solutions of 
the TRILEX scheme above a temperature denoted 
(resp. Tap) for the a;j/z-decoupling (resp. ^-decoupling). 
The evolution of the local vertex and self-energy (resp. 
lattice self-energy and polarizations) is presented in Fig.3 
(resp. Fig. 4) for the points A, B and C of Fig. 2, in the 
xy^-decoupling. At weak coupling (point A), the local 
vertex A reduces to the bare vertex A = 1 at large fre¬ 
quencies, up to numerical noise (Fig. 3a, upper panels). 
The spin polarization (hence the spin susceptibility, see 
Eq.(9)) becomes sharply peaked at the AF wavevector 
Q = (7r,7r) (Fig. 4, upper panels), reflecting the nesting 
features of the Fermi surface. As a result, the self-energy 


































































4 


P‘*(q,w=0) P*P(q,aj=0) ImE(k,iW(,) 



Figure 4: (color online) Momentum-dependence of the self¬ 
energy and polarization (half-filling, t' = Q). A, B and C 
are defined in Fig. 2. Left: ReP‘^*'(q,cu = 0). Middle: 
ReP’**’(q, tu = 0). Right: ImE(k, iwo). 



Figure 5: (color online) From left to right: A(k,tj = 0), 
x(q, = 0) and ImS(k, iwo) in the doped case: U/D = 1.8, 
t' = —0.4t, j3D = 96, 5 = 10%. (lya-decoupling) 

acquires a strong k-dependence at (tt, 0)(Fig. 4), but its 
local part is the same as the DMFT self-energy (Fig. 3b). 
At strong conpling (point C), the vertex becomes simi¬ 
lar to the atomic vertex (Fig. 3a, lower panels). Fnr- 
thermore, the self-energy and polarization are weakly 
momentum-dependent (Fig. 4, lower panels), in agree¬ 
ment with clnster DMFT calculations; the self-energy of 
TRILEX is very close to the DMFT self-energy (Fig. 3b). 
Finally, at intermediate coupling (point B), Ajmp acquires 
frequency structnres which interpolate between A and C 
(Fig. 3a, middle panels), while E is strongly momentnm- 
dependent and its local part departs from the DMFT 
self-energy (Fig. 3b, middle panels). 

The temperature is determined by extrapolat¬ 

ing the inverse static AF snsceptibility (Fig.A.l in Suppl. 
Mat.). It is reduced with respect to the Neel temperature 
computed in DMFT [70] as a result of nonlocal fluctna- 
tions beyond DMFT. Furthermore, ^af- ^ 


consequence of the apparent divergence in the spin sns¬ 
ceptibility at low temperatnres (Fig. A.l), cansed by a 
vanishing denominator of (Eq. 4b), we cannot ob¬ 
tain converged resnlts in the close vicinity of and below 
Taf- Whether we have an actual AF transition or finite 
but very large correlation lengths (as seen e.g in [56]), 
conld be decided by generalizing the present formalism to 
the symmetry-broken phase. Contrary to clnster DMFT, 
the snsceptibilities are not by-products of the calculation, 
but directly enter the self-consistency loop throngh W'^ 
(see Eq.(9)). We thns cannot converge paramagnetic so¬ 
lutions below an AF phase transition. 

Let us now tnrn to the effect of doping. In Fig. 5, 
we present resnlts for t' = —0.4t, pD = 96 and an in¬ 
termediate interaction strength {U = 1.8, close to point 
B). The spectral fnnction displays Fermi arcs (Fig. 5, 
left panel), as observed in experiments [77] and in cluster 
DMFT [35, 37, 38, 42-44, 78] . Let us emphasize that 
this is obtained by solving a single-site qnantnm impu¬ 
rity problem, a far easier task than solving clnster im¬ 
purities. The Fermi arc is a consequence of the large 
static spin susceptibility at the AF wavevector (Fig. 5, 
middle panel), which translates into a large imaginary 
part of the self-energy (Fig. 5, right panel). The corre¬ 
sponding variation of the spectral weight on the Fermi 
snrface is rather mild due to the moderate correlation 
length (^af 2 unit spacings) for these parameters. 

Alternative self-consistency conditions are possible, e.g 
Xk>c = xLp instead of wonld enforce snm 

rules on two-particle quantities that are key to preserv¬ 
ing the Mermin-Wagner theorem in [7, 55] . However, this 
leads to a positive 77 ®p(t) and hence to a severe sign prob¬ 
lem in the quantnm Monte-Carlo at low temperatures. 

In conclusion, we have presented the TRILEX formal¬ 
ism, which encompasses long-range spin fluctuation ef¬ 
fects and Mott physics in a unified way. Like DMFT, 
it can be systematically controlled by extending it to 
clnster schemes that interpolate between the single-site 
approximation studied in this paper, and the exact so- 
Intion of the model. We expect that the convergence of 
the method as a fnnction of the cluster size will strongly 
depend on the decoupling channel and, when done in the 
physically relevant channel, will be faster than cluster 
DMFT methods. Furthermore, because the competition 
between spin fluctuations and Mott physics can be de¬ 
scribed already at the single-site level, the method may 
be a good starting point for correlated mnltiorbital sys¬ 
tems where spin fluctnations play an important role, like 
pnictides snperconductors. 
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Supplementary Materials 

Appendix A: Atomic Vertex 


In the atomic limit (single atomic site), one can compute the three-point vertex exactly by writing its Lehmann 
representation. One gets the following expression [69]: 
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where xlt = § coshf/WO) half-filling, {n^) = 1/2. 


(Al) 


Appendix B: Computation of the three-leg vertex 


Airnp is computed from the fermionic three-point correlation function through the relation: 
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Here, the suffix “c” stands for “connected”, namely: 

xLp,c(*‘^of^) = xLp(*w,iO) - PGimp{i(^){n3^p)Sin 
while is defined as iO) = , i^) ± xltp(*w,i0), with: 


xrmp(*w,iO)= // {Tc^{T)cl{0)n^>iT')) 


x/mp(*0) is the Fourier transform of Xhnp(''’) — {Tn^{T)n'^{0))c- 
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Appendix C: Inverse spin susceptibility: Temperature evolution 

The inverse static AF spin susceptibility (x’^^(Q) = 0)“^, obtained from H/®P(q, iil) by Eq. (9)) decreases linearly 

with temperature, as shown in Fig. A.l. This allows to determine Taf by extrapolation. The resulting phase diagram 
is shown in Fig. 2. 



Figure A.l: (color online) Evolution of the inverse susceptibility with temperature (half-filling, t' = 0). The circles denote 
TRILEX computations, while the crosses on the a;-axis are the extrapolated points. The dashed portion is the extrapolated 
line. Left: 2 -decoupling. Right: a;i/ 2 -decoupling. 
















































